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Abstract 

We investigate cosmological dynamics of multiple tachyon fields with inverse square po- 
tentials. A phase-space analysis of the spatially flat FRW models shows that there exists 
power-law cosmological scaling solutions. We study the stability of the solutions and find 
that the potential-kinetic-scaling solution is a global attractor. However, in the presence 
of a barotropic fluid the solution is an attractor only in one region of the parameter space 
and the tracking solution is an attractor in the other region. We briefly discuss the physical 
consequences of these results. 
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1 Introduction 



Cosmological inflation has become an integral part of the standard model of the universe. 
Apart from being capable of removing the shortcomings of the standard cosmology, the 
paradigm seems to have gained a fairly good amount of support from the recent observa- 
tions on microwave background radiation. On the other hand, there have been difficulties 
in obtaining accelerated expansion from fundamental theories such as string/M-theory. Re- 
cently, Sen 1J has constructed a classical time-dependent solution which describes the decay 
process of an unstable D-brane in the open string theory. During the decay process the 
tachyon field on the brane rolls down toward the minimum of the potential. The tachyon 
field might be responsible for cosmological inflation at the early epochs due to tachyon con- 
densation near the top of the effective potential and could contribute to some new form 
of cosmological dark matter at late times. Several authors have investigated the process of 
rolling of the tachyon in the cosmological background [3 and in the braneworld scenario jl] . 

The stability of tachyon inflation against changes in initial conditions has been studied 
for exponential potentials [H] and for inverse power- law potentials [Sj. Exact solutions for a 
purely tachyon field with an inverse square potential are known [7j, but no solutions exist for 
multiple tachyon fields, so a dynamical analysis may be relevant. It is interesting that the 
inverse square potentials play the same role for tachyon fields as the exponential potentials 
do for standard scalar fields jHj- The potentials allow constructing an autonomous system 
which gives power-law solutions. Recently, it was pointed out that the tachyon model of 
inflation is hard to be consistent with observations unless there is a large number of 
coincident branes [TO] . In this letter, we will consider cosmological dynamics of multiple 
tachyon fields with inverse square potentials. We have assumed that there is no direct 
coupling between the inverse square potentials. The only interaction is gravitational. A 
phase-space analysis of the spatially flat FRW models shows that there exist a cosmological 
scaling solution which is a unique attractor. We find that accelerated expansion of the 
universe can be driven by multiple tachyon fields at lower-Planck energy densities. In this 
model we then introduce a barotropic perfect fluid to the system. We discuss the physical 
consequences of these results. 

2 Multiple Inverse Square Potentials 

We start with more general model with m tachyon fields (pi, in which each has an inverse 
square potential 



Note that there is no direct coupling of the tachyon fields, which influence each other only 
via their effect on the expansion. The evolution equation of each tachyon field for a spatially 
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flat FRW model with Hubble parameter H is 
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where kt = 87rGjv is the gravitational coupling. Defining 2m dimensionless variables 

Xi = 4>i, (4) 
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the evolution equations (J2J) can be written as an autonomous system: 
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where a prime denotes a derivative with respect to the logarithm of the scalar factor, 
iV = In a, and (3i = 4/(3ft 2 Voi). The constraint equation (jH} becomes 
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Critical points correspond to fixed points where x\ = and y[ = 0, and there are self-similar 
solutions with 

This corresponds to an expanding universe with a scale factor a(t) given by a oc t p or a 
contracting universe with a scalar factor given by a oc {—t) p , where 

p = , m 1 ■ (10) 

Setting x'j = and ?/■ = 0, we can get the fixed points K corresponding to kinetic- 
dominated solutions and the fixed points S listed in Table 1 with 
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Table 1: The properties of the critical points in a spatially flat FRW universe containing 
two tachyon fields with inverse square potentials. 



where we using the following definition: 
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The fixed points K behave like non-relativistic matter with p = 2/3. The fixed point S is 
the kinetic-potential-scaling solution with 
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p ~ py (13) 

We find that accelerated expansion of the universe occurs if (3 < 2/y/3. For a single 
tachyon field with an inverse square potential, Vq must be larger than 2/(v / 3k 2 ) to guarantee 
accelerated expansion. However, accelerated expansion of the universe can be driven by 
sufficiently multiple tachyon fields with Vq at lower-Planck energy densities. 

In order to analysis the stability of the critical points, we only consider the cosmologies 
containing two tachyon fields. Using the Friedmann constraint equation (|Sjl. we reduce 
Eqs.flHJ) and (0) to three independent equations 
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The linearization of system (|14j) - (jl6|) about the fixed points K yields three positive char- 
acteristic exponents 

Ai = A2 = 6, A3 = 3, 

which indicate that the kinetic-dominated solutions are always unstable. Substituting linear 
perturbations about the critical point 5* into Eqs. (fT4^) - (flTj) . to first-order in the perturba- 
tions, yields the three negative eigenvalues 
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(2x 2 s + (3 2 ) - y/(2x* + FY - 16^x2(1 + A/AKl - J/fk) ■ 

The critical point is consequently stable so that the corresponding cosmological scaling 
solution is always a global attractor for any (3. For a single tachyon field model (3 = (3i, the 
three eigenvalues reduce to one eigenvalue. 





3 Plus a Barotropic Perfect Fluid 

We now consider multiple tachyon fields with inverse square potentials evolving in a spatially 
flat FRW universe containing a fluid with barotropic perfect equation of state P 7 = (7— l)p 7 , 
where 7 is a constant, < 7 < 2, such as radiation (7 = 4/3) or dust (7 = 1). The evolution 
equation for the barotropic perfect fluid is 

p 7 = -3if(p 7 + P 7 ), (17) 

subject to the Fridemann constraint 
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We define another dimensionless variable z = jffi- The evolution equations (J2J) and (fT7|) 
can then be written as an autonomous system: 
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and the Fridemann constraint equation becomes 
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Critical points correspond to fixed points where x' { = 0, y' = and z' = 0, and there are 
self-similar solutions with 
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Table 2: The properties of the critical points in a spatially flat FRW universe containing 
two tachyon fields with inverse square potentials plus a barotropic perfect fluid. 



This corresponds to an expanding universe with a scale factor a(t) given by a oc t p , where 



The system fll9)) - (j21|) has a fixed point K corresponding to kinetic-dominated solution, 
a fixed point S which is a kinetic-potential-scaling solution, a fixed point F which is a 
fluid-dominated solution, and a fixed point T which is a tracking solution listed in Table 
2. In order to analysis the stability of the critical points, we still consider the cosmologies 
containing two tachyon fields plus a barotropic perfect fluid. 

K: Xi = ±1, y = 0, z = 0. These kinetic-dominated solutions always exist for all Pi 
and 7, which behave like non-relativistic matter with a oc t 2 ^ 3 irrespective of the nature of 
the potentials. The linearization of system ([19 )1 - ()21 )) about these fixed points yields four 
eigenvalues 

Ai = A 2 = 6, A 3 = A 4 = 37, 
which indicate that the solutions are always unstable. 

S: Xi — x s , y — x 2 s /Pi, z — 0. The potential-kinetic-scaling solution exists for all Pi 
and 7. The power-law exponent depends on the parameter P of the potentials. The 
linearization of system -(|2*T |) about the fixed point yields four eigenvalues 



which indicate that the solutions are unstable for 7 < j3(y//3 2 + 4 — (3)/2 and stable for 
7 > P(y/p* + 4-/3)/2 from Eq.flUJ). 

F: Xi — 0, yi — 0, z — 1. The fluid-dominated solution exists for all /3j and 7, 
corresponding to a power- law solution withp = 2/37. The linearization of system ()19j) - (j2*T]) 
about the fixed point yields four eigenvalues 

Ai = A 2 = -3, A 3 = A 4 = 37, 

which indicate that the solution is unstable for 7 7^ and stable only for 7 = 0. 

T: Xi = ^7, yi = j//3i, z = 1 — pj[ ■ The tracking solution exists for 7 < j3(y/(3 2 + 4 — 
f3)/2. The solution displays a tracking behavior according to the definition in [H]. The 
power-law exponent, p = 2/3-y, is identical to that of the fluid-dominated solution, depends 
only on the barotropic index 7 and is independent of the parameters (3i of the potentials. 
The linearization of system (fT9"|) - ()21)) about the fixed point yields four eigenvalues 
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which indicate that the solution is always stable when this point exists for 7 < (3{^/ (3 2 + 4 — 
(5)12. 

The different regions in the (/?, 7) parameter space lead to different qualitative evolution 
in Figure 1. For 7 > /3(y/ j3 2 + 4 — /3)/2, if, S and F exist. Point S is the stable late- 
time attractor. Hence generic solutions begin at the two kinetic-dominated solutions or 
at the fluid-dominated solution and approach the kinetic-potential-scaling solution at late 
times. For 7 < (3(y/ '/3 2 + 4 — /3)/2, all critical points exist. Point T is the stable late- 
time attractor. Hence generic solutions start at the two kinetic-dominated solution, at the 
kinetic-potential-scaling solution or at the fluid- dominated solution and approach the stable 
fluid-kinetic-potential-scaling solution. 

4 Conclusions and Discussions 

We have presented a phase-space analysis of the evolution for a spatially flat FRW uni- 
verse containing m tachyon fields with inverse square potentials. We find that there exist 
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Figure 1: Stability regions of the 7) parameter space. In the regions I and II, the 
tachyon kinetic-potential scaling solution is the stable late-time attractor. In the regions 
III and IV, the fluid-dominated solution is the stable late-time attractor. The universe 
accelerates in the regions I and III, while the universe decelerates in the regions II and IV. 



cosmological scaling solutions that corresponds to an expanding universe with a oc t p . For 
a single tachyon field with an inverse square potential, the universe could accelerate only 
at nearly Planckian energy densities. However, accelerated expansion of the universe can 
be driven by sufficiently multiple tachyon fields even at lower-Planck energy densities. The 
reason for this behavior is that while each field experiences the 'downhill' force from its 
own potential, it feels the friction from all the tachyon fields via their contribution to the 
expansion \1'2\ . In order to analysis the stability of the critical points, we only consider the 
cosmologies containing two tachyon fields. We find that the critical point is always stable 
so that the scaling solution is a global attractor irrespective of the form of the potentials. 
This implies that the velocity of each tachyon field tends to be equal and constant via their 
effect on the expansion. 

Then we have extended the phase-space analysis of the evolution to a realistic universe 
model with a barotropic perfect fluid plus m tachyon fields with inverse square potentials. 
We have shown that the energy density of the tachyon dominates at late times for 7 > 
f3(Vf3' 2 + 4 - (3)/2, In constraint, for 7 < f3(y//3 2 + 4 - (3)/2, the barotropic fluid does not 
dominate completely and the contribution of tachyon energy density to the total one is not 
negligible. 

We emphasize that we have assumed that there is no direct coupling between these 
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inverse square potentials. In a system with m non-coincident but parallel non-BPS D3- 
branes [IB], there are two kinds of open strings. One of them starts from and ends on 
the same brane; The other starts from a given brane, then ends on a different brane. If 
the distance between two branes are much larger than the string scale, one can ignore the 
second kind of open strings, leaving a tachyon on the world volume for every brane. Thus, 
we have m tachyons without interaction and the action is simply the sum of m single- 
tachyon actions. It is worth studying further the the cosmological dynamics of multiple 
tachyon fields with interactions. 
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